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We compute the amplitude of the non-Gaussianities in inflationary models with multi- 
ple, uncoupled scalar fields. This calculation thus applies to all models of assisted infla- 
tion, including N-flation, where inflation is driven by multiple axion fields arising from shift 
symmetries in a flux stabilized string vacuum. The non-Gaussianities are associated with 
nonlinear evolution of the field (and density) perturbations, characterized by the parameter 
Jnl- We derive a general expression for the nonlinear parameter, incorporating the evolu- 
tion of perturbations after horizon-crossing. This is valid for arbitrary separable potentials 
during slow roll. To develop an intuitive understanding of this system and to demonstrate 
the applicability of the formalism we examine several cases with quadratic potentials: two- 
field models with a wide range of mass ratios, and a general A/"-field model with a narrow 
mass spectrum. We uncover that /at l is suppressed as the number of e- foldings grows, and 
that this suppression is increased in models with a broad spectrum of masses. On the other 
hand, we find no enhancement to fj^L that increases with the number of fields. We thus 
conclude that the production of a large non-Gaussian signal in multi-field models of inflation 
is very unlikely as long as fields are slowly rolling and potentials are of simple, quadratic 
form. Finally, we compute a spectrum for the scalar spectral index that incorporates the 
nonlinear corrections to the fields' evolution. 
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INTRODUCTION 



d models of inflation have been considered ever since the introduction of hybrid infla- 



In these models, only one field typically evolves during inflation. The role of the 
second field is to add a potentially tachyonic direction to the potential which ends inflation by 
creating an instability in a direction orthogonal to the classical inflationary trajectory. Conversely, 
assisted inflation [4J relies on J\f uncoupled fields. The fields are each unable to generate a workable 
period of inflation on their own, but can evolve coherently to provide a cosmologically acceptable 
inflationary epoch. 

The analysis here applies to any model where inflation is driven by multiple scalar fields that do 
not interact directly with one another. This restriction is less onerous than it might appear, since 
assisted inflation typically only works to the model builder's advantage when the individual fields 
are very weakly coupled to one another [5J . Our treatment of this problem is motivated by N-flation 



[fi[Q] 1 - This is a proposal to implement assisted inflation in the context of string theory, and the M 
fields arise from axion potentials associated with shift symmetries in the compact manifold. From 
the model-building perspective, this provides a mechanism for generating inflation without invoking 
field(s) whose VEVs become trans-Planckian at some point during their evolution - a point at which 
perturbative descriptions of any string or supergravity derived model will generically break down. 
Moreover, the shift symmetries of the axions suppress couplings between the fields, a prerequisite 
for successful assisted inflation. It is not yet known whether the microscopic physics of string theory 
will permit N-flation to naturally arise in a realistic scenario, but it certainly represents a novel 
and interesting approach to string inflation. Originally, the axion fields that drive N-flation were 
assumed to have identical masses [6]. Moreover, the sinusoidal axion potentials were approximated 
by Taylor expansions about their minima, retaining only the quadratic term. Subsequently, Easther 



a, 



and McAllister showed that the mass spectrum could be derived via random matrix arguments 
avoiding the intractable calculation required by a direct assault on the problem. The resulting 
distribution of masses conforms to the Marchenko-Pastur distribution, and is controlled by a single 
parameter (3 - the ratio between the number of axions M and the total dimension of the moduli 
space. 

In the limiting case where all the masses and initial field values are identical 2 the spectrum of 
scalar and tensor perturbations from N-flation is identical to that generated by a single massive 
field. Since multi-field models of inflation generically have a richer phenomenology than single-field 
examples, we might hope to break this degeneracy by looking for specific signatures of multi-field 
evolution within N-flation. In particular, the non-Gaussianities in single field models are typically 
ll,Q,[l3] , but they can be larger in general multi-field models, so computing them is an 



tmy 



9 



10, 



11 



important step towards developing a full understanding of N-flation 3 . Several studies have argued 
that the non-Gaussianities from assisted inflation are small [if], Q, [l^l, but these calculations 

n 

contain significant simplifying assumptions. For example, only two fields are considered in [16Q . 



and 



13, 



181 ] use the horizon crossing approximation. As a result, a general expression for the 
non-Gaussianities, described by /nl, is lacking. To be more precise, terms are usually neglected 
that incorporate the evolution of perturbations once they cross the horizon and the effects of 



1 See also [g] for the derivation of a closely related assisted inflation model based on multiple branes, to which our 
study also applies. 

2 When assisted inflation is driven by m 2 (f> 2 potentials the inflationary dynamics have some "memory" of the initial 
field values, whereas other widely studied implementations of assisted inflation have an attractor solution in field 
space. 

3 An additional source of Non-Gaussianity is the era of reheating, which we do not consider in the following; see eg. 
[13 . UBI ] and references therein for relevant literature. 
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isocurvature modes (both of which are always possible in multi-field models) - see for example 

18]. 

We derive this desired general expression within the <5 TV-formalism, first proposed by Starobinsky 
fig| and extended by Sasaki and Stewart 20] among others [if], Q]. The only assumptions we make 
are that the potential is separable and that the slow roll approximation is valid. Since the general 
expression contains a part that is not immediately slow roll suppressed, we will focus our efforts 
on this term. To illustrate the application of the formalism and to build up intuition we consider 
several specific cases with quadratic potentials: first, an exactly solvable two-field model where the 
ratio of the squared masses is two; second, generic two-field models where the ratio of the squared 
masses is less than five; third, the generic multi-field case with a narrow mass spectrum. 

In the two-field cases, we find an unexpected suppression of Jml by the volume expansion rate, 
expressed in terms of the number of e-foldings. The exponent of the rate is given by twice the 
square of the ratio of the heavier mass to the lighter one. Based on this result, we focus our 
attention on narrow mass spectra, for which we compute the nonlinear parameter up to second 
order in the width of the spectrum (properly defined in Sec. IHICj) . As expected, we encounter a 
suppression factor that scales as the square of the number of e-foldings. However, we do not find 
any enhancement that scales with the total number of fields. Possible exceptions to this argument 
include models where one or more fields violate slow roll, or where there are significant couplings 
between the fields. We intend to address these issues in future work. However, at present we find 
no evidence that the non-Gaussianities generated by assisted inflation modes - including N-flation 
- are enhanced relative to those of their single field analogs. 

The article is structured as follows: first, we derive the general expression for /tvl in Sec. Ill Bl 
followed by a discussion where we reduce this general result to recover the specific cases already 
encountered in the literature. Thereafter, we discuss three specific cases, two-field models in 
Sec. IIII Al and HTl Bl and the multi-field model in Sec. lIIIC| again followed by a discussion. Last but 
not least, we derive the general expression of the scalar spectral index incorporation the evolution 
of perturbations once they cross the horizon in Sec. but we postpone specific case studies to a 
follow-up publication. We conclude in Sec. [V] 

II. NON-GAUSSIANITIES IN MULTI-FIELD INFLATION 



We are interested in evaluating non-Gaussianities in models of inflation with multiple scalar 
fields, in the hope of finding possible experimental signals that will break the degeneracy between 



5 



multi-field models and common single field models. 4 We shall use the <5iV-formahsm to compute 
/jVX> which characterizes non-Gaussianities. This formalism was prop osed by Starobinsky [lj]] 



16 



13, 



2l\ . In this approach 



and further developed by Sasaki and Stewart in [20|], and others in 
one relates the perturbation of the volume expansion rate 6N to the curvature perturbation £ 
on large scales, which is possible if the initial hypersurface is flat and the final one is a uniform 
density hypersurface [20]. Note that £ is conserved on large scales in simple models, even beyond 



linear order 



25 



261 ] 5 . Given this relationship between the curvature perturbation and the volume 



expansion rate, one can evaluate the momentum independent pieces of non-linear parameters, 
which are related to higher order correlation functions, in terms of the change in N during the 
evolution of the Universe, see e.g. j^ . 

Our treatment will parallel that of Vernizzi and Wands [16], who computed the general ex- 
pression for /jvl, but restricted themselves to two-field models when they came to compute the 
non-Gaussianities. 

To begin, consider the action for N scalar fields, 



m 2 r (\ N \ 



(1) 



which we assume to be responsible for driving an inflationary phase. The unperturbed volume 
expansion rate from an initial flat hypersurface at t* to a final uniform density hypersurface at t c 
is given by 



Hdt. 



(2) 



where H is the Hubble parameter. The nonlinear parameter /jvz, can be related to the derivatives 
of the expansion rate N with respect to the field values ipj(t*) = tp*. This computation, starting 



from the three point correlation function, was performed in 



16 



171 ] resulting in 



5 fNL 



Z =1 NiNjN^ 



where the short hand notation 



Ni 



dN 



(3) 
(4) 

(5) 



4 See [171 . l2ll 1 2 21 . |23| for a sample of the recent literature on multiple field models and e.g. [9|, ll(j . |24| for a discussion 
of single field models. 

5 The separate Universe formalism put forward by Rigopoulos and Shellard in e.g. [2^] is equivalent to the 8N- 
formalism. 
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d 2 N 

N * S d^* ' (6) 
was used (we refer the reader to [la, for details). The first term in ([3]) is small. On geo- 

nn 

metrical grounds, we know < / < 5/6 [3, llq |. while r is the usual tensonscalar ratio. 6 The 
observational upper limit on this quantity depends on the priors used in the fitting process, but we 
can reliably conclude that r/16 < 0.1 28(. Observationally, it is very unlikely we will ever detect 
non-Gaussianities unless Jnl > 1- Henceforth, we focus on the second term in ([3|), to determine 
under what conditions non-Gaussianities could become large. Currently, the best observational 
bound on f NL is drawn from the WMAP3 data Q: -54 < f NL < 114. Recently, Q Kim and 
Liddle derived the estimate (6/5)/]$ < r/16 by ignoring the evolution of perturbations once they 
cross the horizon, and constraining the potential W to be the sum of monomials in (p^. Conversely, 



161 ] is restricted to the two-field case but includes effects from evolution that occurs after horizon 



crossing, and arrives at a similar conclusion. Since sea' 
horizon crossing in multi-field models, we extend 16], 



ar perturbations need not freeze out after 
18 1 to compute /]$ for N fields without 



assuming any freeze-out. Our principal assumption is that the slow roll approximation is valid for 
all fields - in practice, individual fields can cease to be critically damped well before inflation comes 
to an end, depending on their masses and initial values. Since we retain any possible super-horizon 
evolution, our analysis follows a similar path to [l6]. 

A. Slow Roll Dynamics 

In the following, we will set the reduced Planck mass m p = (87rG) -1 / 2 = 1 for notational 
simplicity. Since we are interested in assisted inflation, we can ignore cross coupling terms between 
the scalar fields. That is, we assume 

N 

W{^^ 2 ,...) = Y,Vi^i), (7) 

i=l 

but we keep the form of the potentials Vi(ipi) general. The equations of motion for the fields are 

(pi + 3H<pi + V( = , (8) 



6 The quantity / incorporates the dependence of the three point correlation function on the shape of the momentum 
triangle; the maximal value results for an equilateral triangle and the minimal one if two sides are much longer 
than the third [jj. 



7 



where we set V- = W% = dVi/dcpi. The Friedman equations read 



H 2 



i=l 

Slow roll inflation occurs if the slow roll parameters, defined as 7 

1 V/ 2 VI' 

£i = 

are small (e; < 1, jjj < 1), and 



H = -Y.\$- (10) 



AT 

e = J^£i«l. (12) 

j=i 

In this case the dynamics is governed by 

3H<pi « -F/, (13) 
3# 2 » W. (14) 

For simplicity, we assume V- > from here on, so that we can replace the derivatives by the slow 
roll parameters, V( = W \f2e- L . 

During slow roll inflation we can write the number of e-foldings ([2]) as 

rc AT y 

N(t c ,U) = - / J2^difi. (15) 

i=i » 

Furthermore, there are — 1 integrals of motion, for example the set 



dtpi f dip 



c — i? + ifr' (16) 

for i = 1 ... M — 1. These Cj can be used to discriminate between different trajectories in field space 
and they will become quite handy in the next subsection when we evaluate iVj and N^j. 



B. Non-Gaussianities 

We will now evaluate the derivatives of the volume expansion rate with respect to the fields, 
which are needed in order to evaluate the nonlinear parameter. First, write down the total differ- 
ential of N by using its definition (|15p . 

dN =Yl 



3=1 



V / ^ dip* \ V 



dip* . (17) 



7 We are using the potential slow roll formalism here. For a more general treatment of the slow roll expansion with 



multiple fields, see [29(. 
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Furthermore, using the integrals of motion Cj from (|16p . we can relate d<p? and dip* via 

*5 = Eai(Es>i). < 18 > 



where 

= y=^(^-i"^)- ( 19 ) 

Next, we would like to make use of the M — 1 integrals of motion to eliminate d(pj/dC{ in favor 
of d(fi/dCi, which can then be used in the condition p = const at t c , the time at which we want 
to evaluate the non-Gaussianities. To accomplish this, consider 

i-l 

C, = ( 2 °) 
Differentiating this with respect to Ck yields 



l j i 



which can be solved to give 



da = _dtf_^_dtfj_ 

dc k dc k (y{)c dc k (vf) c ' { ) 



where we introduced 



[ 1, if Jfc < i - 1 
&ki=< ~ (24) 

[ 0, if k > i - 1 . 

If we plug this into the derivative (with respect to C{) of the p = const condition, 



after some algebra we arrive at 



S-^'fef)/ 9 ^ (26) 



( E e , c - e ^ c I ■ ( 27 ) 

U=fc+1 



9 



where we used the definition of the slow roll parameters (jlip in the last step. Using (|27p and (fl~9 
in (|18j) we end up with 



(28) 



We are now ready to compute the derivatives of the expansion rate from (|17j) , which reduce to 

(29) 



aiv _ i v* + z= 



where we introduced 



i=l 



W c ^r - V k c . 



(30) 
(31) 



After some algebra, the second derivative becomes 



with 



d 2 N 
MM 



2e* W* 



dZf 



'—W*A lk 



+ 



1-S 



(32) 



(33) 



(34) 



Note that Aik is a symmetric matrix. Given that, we can now write down the general expression 
for ffil as 



6 f (4) _ ^=l eg y 1 ^ 2e* ; + 2^fe,i=l e ^ ^ 
5 / iVL ~~ z 



2 N 2 



where is given by 



Uk 



AT4 



(35) 



(36) 
(37) 



with AV fc = V* - v§ > 0. 



Equation (|35p is the desired nonlinear parameter characterizing non-Gaussianities and our first 
main result. The factor Z c incorporates the evolution of perturbations after they cross the horizon 
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until the time t c - it is this factor that was neglected in 18]. Before we discuss (|35p in the 
next subsection, it is worthwhile to reiterate our assumptions. First of all, we assume slow roll 
inflation throughout. This is a tricky assumption in inflationary multi-field models, since the 
heavier fields can leave slow roll while inflation still continues, see e.g. Qj. In this case, computing 
the contribution of these fields onto /jvx would require us to go beyond the formalism used here. 
Next, we assumed a separable potential, neglecting cross coupling terms between the fields, but 
otherwise we kept the form of the potential general. Neglecting cross couplings is reasonable 
when examining assisted inflation models like N-flation, since these models require that fields do 
not interact to any significant degree. Last but not least, some minor technical assumptions are 
introduced via our use of the (5iV-formalism, such as the reduction of four-point functions to a 
product of two-point functions by means of Wick's theorem - we refer the reader to 
for the derivation of the e)./V-formalism. 

C. Discussion 

Our main result is the nonlinear parameter given by (|35p together with (|34p . In the case of two 
fields, this reduces to the case examined by Vernizzi and Wands [16]. To be specific, with M = 2 
equation (j34|) becomes 

•An = A 22 = - A12 = -Ai\ = A , (38) 



where we introduced 



W} e. 

As a consequence, (|35j) becomes 



A^^Mf^mim) . (39 ) 



2 / „* \ „2 



^l-^i)+^(l-lH + (^-f ) A 

5 JJVi - " /..-' „2\2 



6,(4, _ 2 ,-A ■ w _ (40) 



2 



which is identical to (68) of 16] 



a 1 a 

* T F * 



On the other hand, if Z£ ~ we recover the simple expression of jl8| for M scalar fields, which 
resulted in a small contribution of f$ L to the total nonlinear parameter. 

Consequently, a large non Gaussian signal could only arise if Z? becomes large. A closer look 
at (|35p reveals that only the last term in the numerator is not immediately slow roll suppressed, 
that is the term proportional to the ^4-matrix elements defined in (|34p . Therefore, let us look at 
this matrix and ask under what conditions its contribution could become large. 
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First, note that the sum over a column or a row of the ,4-matrix vanishes, since 



Af 



1=1 



0, 



M 



E%(f-^-?)E(?-% 

j=i 1=1 



(41) 



because E^ E\ = e so that the last sum is identical to zero. Based on this fact one can immediately 
see that the contribution to the ^4-matrix vanishes if we deal with J\f identical fields (as expected). 
But even in the generic case, we should expect some cancellation. 

Next, let us go back to the definition of Uk in (f3TJ) : if AVk <C W c e c k /e c , we can approximate 

W c el 



ui. ~ . 

W* e c 

Using this in (134p . we can write the last term in (135D as 

2 



(42) 



UkU l A„ 
l^k,l=\ els*^ 

l^k=l el _ 



£2 



£1 



1 + 



+ 2er% 



V2 



(43) 



where we introduced 



c\2 



f (£|) 
k=i k 



Af 



c\3 



£2 



y (4) 
^ (4) 2 



(44) 



k=l 



and 



Af 

E 

fc=i 



AT 



Af 

E 



c\3 



(4) 



c( c *\2 



k=l ' ^* fc=l 6 ^ 

As an application, let us have a look at the "horizon crossing limit" t c — * i*, which corresponds 
to purely adiabatic perturbations. In this limit e~i = e* = e c and fji = E^e^r^/e*, so that (fl3l) 
vanishes identically. One can see this from (|35j) . since the prefactor in front of the ^4-matrix 
elements becomes independent of k, so that a sum over the matrix elements is identically zero 
due to (|41|) . The second term in (|35|) is sensitive to the evolution of modes after horizon crossing. 
This evolution and thus the non-Gaussianity is closely correlated with the presence of isocurvature 
perturbations. These modes generically occur in multi-field inflationary models whenever light 
degrees of freedom transverse to the adiabatic direction are present, but it is difficult to transfer 
them to the adiabatic mode during slow roll inflation: in order to transfer isocurvature modes 
effectively, the trajectory in field space should be sharply curved, which might occur when a field 



Vk 



(45) 
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leaves slow roll. However, this process can not be described properly with the formalism at hand, 
since the slow roll conditions are necessarily violated. In addition, we can generate isocurvature 
modes at the end of inflation if the different fields decay into different particle species. 

We now look at another interesting case for which our formalism is valid, and which might still 
permit non-Gaussianities, even though we do not expect them to be large. Assume the time t c lies 
towards the end of inflation, just before slow roll ends. Even if ft\ , £2 > e c *C 1 we could still get a 
non-Gaussianity of order one if the r\i are reasonably large. First of all, it seems possible to cancel 
the suppression due to £%/e c in fji, which is of order 0{\/N): since we have N summands, we can 
naively expect an enhancement of order 0(N) that cancels out the previous suppression. Next, if 
we are near the end of inflation, we have e£ S> £ Zi so that the ratio becomes large. However, 

this enhancement factor in fji and 772 gets compensated by the prefactor. Thus, we might naively 
expect fffi L = 0(770) towards the end of inflation, which can be of order one. 

However, we have to be more careful here: we already saw in the case of M identical fields that 
cancellations occur in the sum, due to the symmetries of the ^4-matrix. In the next section we will 
examine a few specific models in order to develop intuition as to whether non-Gaussianities can 
become large. It will turn out that the naive argument above is misleading, since the sum over 
the ^4-matrix generically leads to terms that are slow roll suppressed. Hence, it seems very hard 
to produce a significant non-Gaussian signal in assisted inflation models which are well-described 
by slow-roll. 



III. CASE STUDIES 



In section III Bl we derived the general expression for the nonlinear parameter, (|35p . The only 
term which was not obviously small during slow roll is 

Fs ?±tiM±. (46) 
(^4) 

Consider this expression for a few specific models: first, an "exact" solvable two-field model with 
ml = 2m^; second, approximate two-field models with m\ = amf and 1 < a < 5; third, AA-field 
models with m\ / mf = /x? = 1 — <5j where 5i <C 1 . 

The first step consists always of evaluating the field values at and t c , which we denote by ip* 
and To do so, we need 2J\f conditions, which are given by: N — 1 dynamical relations between 



In N-flation, this happens frequently during inflation, since the heavier a field is the earlier it will leave slow roll 

a. 
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the fields, N — 1 initial conditions, one condition from the requirement that be N e-foldings 
before t c , and one condition from the requirement that slow roll ends for at least one field at t c . 
Thereafter, we need to evaluate the slow roll parameters at t* and t c , which in turn enable us to 
compute F. 

Consequently, let us first look at the dynamics of multi-field inflationary models with 

JV 

W = Y,Vi (47) 
i=i 

= lEm^l (48) 
z »=l 

where we ordered the fields such that rrn > rrij Hi > j. The field equations (fTB"|) during slow roll 
yield 

where we also used the Friedman equation (|14p. Dividing equations such that W drops out and 
integrating yields the J\f — 1 conditions 

5- (If- 

It should be noted that this is not an attractor: initial conditions for the fields influence the 
dynamics at all times. As a consequence, a thorough study of initial conditions should be performed 
whenever a multi-field model with separable quadratic potentials is proposed as a serious scenario 
of the early universe. However, we are not interested in testing a particular model right now, but 
in developing intuition by studying a few concrete examples. Naturally these models are chosen 
such that they are easy to treat. Therefore, we restrict ourselves in most cases to equal energy 
initial conditions, that is we impose V* = Vf, which will simplify matters considerably. Of course 
this results in the N — 1 initial conditions 

<Pi = IHVx ■ (51) 

Next, consider the number of e-foldings defined in (|15p : this expression integrates to 

AT 

AN = £ [(^*) 2 - (tf) 2 ] ■ (52) 

i=l 

The last missing equation is due to the fact that slow roll ends for one of the fields at t c , that is 
either one of the slow roll parameters rji, £j defined in (|11|) or e becomes of order one. We can 
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easily compute the ratio of the r^'s, yielding rji > rjj Hi > j. Similarly, the ratio of rji to e, becomes 
W/Vi, which is larger than one so that rji > £j. Consequently, the field with the largest mass, that 
is the AA'th field, will leave slow roll first when 

^ = 1, (53) 

given that e c < 1, which will be satisfied in all cases we study. 

Equations (f50|) - (|53|) are the 2J\f conditions which are needed in order to evaluate all fields at t* 
and t c . For any concrete models we can solve the above conditions and evaluate F. For simplicity, 
we suppress the superscript c in the following. 



A. Two fields: mf/m| = 1/2 



We start with this case, since it is possible to solve (150p -([53j) without any approximations. 
Equation ([53]) immediately yields 

tpl = 2(2 - yl) , (54) 

and (fSTj) becomes 

rf 2 = 2<p? . (55) 



Using these in (|52j) yields 



3 

which in turn can be plugged into (|50p with the result 

(2 ~ V?j) 2 _ Q , 
4(iV + l)-^ 

This quadratic equation in Lp\ can be solved to 



9 ^:;, w - 2 = 3^ . (s?) 



(4 = 2 + + + N* (58) 

where we expanded in terms of 1/N in the last step - we need to keep terms up to order iV~ 3 , to 
recapture properly the leading order contribution to F. The remaining fields become 

vf = \{N + l)-^ + ^-^ + 0(N-*). (62) 
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It is now straightforward but somewhat tedious to evaluate all slow roll parameters as well as u\ 
and v,2- In the end, we arrive at 

F = -Tim + ° {N ~ 5) - (63) 

The important feature is the obvious slow roll suppression F ~ N~ 2 ' 2 . This suppression is due 
to the cancellations that occur within the summation over the *4-matrix. As a result, A defined in 
(|39p becomes proportional to iV -4 , so that even though the prefactor in (|40p is of order one, the 
final result is heavily slow roll suppressed. 

B. Two fields: mf/m| = 1/a 

This case can still be solved approximately and we will be able to see how F, and with that the 
non-Gaussianity, depends on the ratio of the masses. Taking analogous steps to the discussion in 
the previous section we arrive at 

pi = a(2 - ifl) , (64) 

<Px = ®vf , (65) 
V? = 4iY + 2a + f (1 - Q) . (66) 

Using these, and writing x = ip\ gives 

{2-x) a (l + a) a 



x(l + a) , (67) 



(4iV + 2a + x(l-a)) Q " 1 
which could be solved numerically. However, we refrain from doing so since we are primarily 
interested in the analytic form of the slow roll suppression of F, which can be computed if we make 
some minor approximations. First, assume that (fi <C (p*- with the consequence that (|51|) and (|52|) 
immediately yield 

*2 _ 4AT 
1 + a 

(70) 



Equation (i53|) still leads to 



p\ = a(2 - <p\) , (71) 

(72) 
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and (f50|) becomes 

a-1 



{2-x) a ( 4N 



a + 1 / 

If we now Taylor expand the left hand side for i <C 1, we arrive at 



(73) 



1_ H + (JO-)" 1 , (74) 

where we truncated the expansion such that a quadratic equation for x results. Note that the 
above expression yields the exact result in case of a = 2. Solving this equation results in 

a—l / r . »r 1 a— I 



2 a ~ 1 a + 



4JV 

Q + l 



A a - 1 a(2 - a) + 



4iV 



a+l 



2 a a + 



4iV 

Q + l 



2a-2 



W — V 2^a(a - 1) L — ■ (75) 

As in the previous section, we can now expand the field values to any desired order in 1/N, compute 
the slow roll parameters, evaluate u± and 112 and in the end compute the non-Gaussianity due to 
F. Doing so we arrive at 

F ~^> ( 76 ) 

for a = 2,3,4, and F = for a = 1. The suppression is again due to cancellations within the 
.A-matrix. Even though the proportionality factor increases with a, it is easily dominated by the 
additional suppression due to N~ 2a . One should note that our assumption ipi <C y*. breaks down 
for a ~ 5, since ip\ ~ v^ 2 /^ already, if we use N = 60. 

We conclude that non-Gaussianities become more increasingly suppressed as the mass difference 
increases. If the masses are similar (a ~ 1), we should still have a suppression of 1/N 2 . Of course 
the non-Gaussianity due to the ^4-matrix vanishes identically if a = 1, that is if the masses are 
equal. Hence, we predict the largest non-Gaussianity to occur for nearly equal masses. To be 
specific, if we write mf/m 2 , = 1 — 5 with 5 < 1, we expect F ~ 5@ /N 2 with some exponent (3. As 
soon as 5 becomes of order one, the additional slow roll suppression kicks in and suppresses the 
signal. 

One assumption we made was the equal energy initial condition (I65p . and one might wonder 
how a different choice effects F: it turns out that only the numerical prefactor in (|76|) changes if we 
choose different initial values for the fields (e.g. for a = 2 it changes from —9/16 to —4 if ipl = <p^ 
is used), but the important suppression factor N~ 2a remains unaltered. As a consequence, we will 
restrain ourselves to the equal energy initial condition in the next section, where we will be able 
to derive the exponent (3 and the proportionality factor in general for multi-field models. 
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C. TV fields: ml/mf = 1 - 5 l 

Let us now consider the case of N fields and focus on a narrow mass spectrum, since a broad 
distribution is expected to result in heavily suppressed non-Gaussianities, based on the analysis in 
the previous section. Another reason to focus on narrow mass spectra is that broad distributions 
are not well suited for assisted inflation, since the heavy fields would roll quickly to their minimum 
without contributing much to inflation. 

Therefore consider 

4 = ( 77 ) 

mf 

with <5j <C 1 and Si > Sj if z > j, that is we order the fields according to their masses. We 
characterize the width of the mass distribution by 

1 M 

i=l 

and 

jfhi p») 



s 2 



Our first task is again to find the field values at and t c via (|5(J| ) - (|53| ). What simplifies matters 
is that we can expand in terms of the <5j. However, as we shall see later on, terms linear in 5i cancel 
out exactly in the expression for F. Therefore we keep all terms up to 5f right from the start. 

By using (i50j) and (f5TI) in (f53l) we get 

and (|52p turns into 

4iV = - 5,) - ^ 2 f> - 5,) (4) " ■ (81) 

i=i i=i V^i / 

without expanding anything so far. Next, noting that 

(f ) ^ - ^ (i + <* + « ) ) + «5 '» 2 ($ ) + °«>) < 82 » 



and using the definitions of (5 and 5 2 , we arrive at 

■¥.<p? » (27Y + 1) +^ + 5^ + 1 -pf 
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+fir ~ <fi y + & ^ (^) ) + ss M + ^ 2 (2iv + 1) , (83) 



^2 ^ 1 + ^_ Jln ^ 

- (5 2 + ^) In + (V - ^ In 2 (^) , (84) 

where we kept all terms up to second order. To solve these coupled equations, we iteratively insert 
them into each other until we have all terms up to second order. The solution reads 

—<p? « 2N + 1 + cW + 2N5 + (5 2 - <5 2 ) 7 + 5 2 M + 2iV,5 2 , (85) 

N 2 ic ff 2iV ff/ff r n 
y^ 2 « i+^ + 7 «j+-^_ -<f(<J-<W) 

+J 2 - + (5 2 -5 2 + ^) 7 + ^^ 7 2 , (86) 

where we introduced 

7 = m(2iV + 1) , (87) 

which is of order one. All other fields follow directly from (|50p and (|51|) . which have to be properly 
expanded too. 

Now that we have the field values, we are ready to compute the slow roll parameters. The 



easiest one is rji = mf/m'ir, which becomes 



« l-6tf + 6i-6tf6i + %. (89) 

Next, we would like to compute e = since the ratio of rji and e appears in the .A-matrix. 

Noting that £» = r}iVi/W as well as TV = mjy-, and plugging in (^ 2 by using ([50]) and (f5"T|) with y? 2 
and f^ 2 from above we arrive after some tedious but straightforward algebra at 

e^l-5 N + 5-55M+~5 2 (i + l)+5 2 1 . (90) 

Note that e < 1 since <5tv > 5, so that 77^ is indeed the largest slow roll parameter. The ratio we 
are interested in becomes 

1 - ^ « 5 - ^ - 5 2 + 5i5 + <5 2 ( 7 + 1) - ,5 2 (1 - 7) . (91) 

At this point we should step back for a second and have a look at the expression for F: first, 
note that X^(<5 — Si) = 0. Because of this and since the leading order contribution of 1 — r/j/e is 
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already first order in delta, we know that F is identical to zero up to first order in delta. This is 
the reason why we computed 1 — rji/e up to second order. Next, since 1 — r/j/e has no zeroth order 
contribution, we only need to evaluate the remaining constituents of F up to first order in delta, 
which simplifies matters quite a bit. 

Doing so leads after a lot more straightforward but tedious algebra to 



(2JV + l) 2 V M(2N + 1) 
25 2 ln(2JV + l) 



1 

AT 2 



1 



2N + 1 



- — — 

M + M 2 



(92) 



(2iV + l) 2 

This is our final result, but before we discuss its implications, it might be instructive to see how at 
least the leading order contribution in the number of fields comes about. By looking at the general 
expression for F in (|46p as well as the definition of the .A-matrix in (|34p . one can convince oneself 
that 

-2 

W 2 I u% \ 




(93) 



3=1 



where \o(8P) means that the adjacent quantity has to be expanded up to order in all deltas. To 
evaluate this expression we use 1 — rjj/s from (f9Tj) . e from ([90]) and 



£ k 
u k 



in 



w 2 



e>(<5°) 



^(l-**(7-l)-<W + 7*) , 

2N+1, 

2N + 1 



1 



(94) 
(95) 

(96) 

(97) 



o(so) (2N + l) 2 ' 
which leads to the O(J\f ) contribution to F from above. 

There are two interesting features to our result ([92]): first, the expression is again suppressed by 
the number of e- foldings (oc N~ 2 ), just as expected from our experience with the two-field cases. 
Next, the leading order contribution is not enhanced by the number of fields, as one might naively 
expect, but is of order 0(5 2 ). As a consequence, the contribution to the non-Gaussianity due to 
the ^4-matrix is negligible, even for multiple fields. 
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We restricted ourselves to a narrow mass spectrum here - however, for a broader mass spectrum 
we expect an even stronger suppression by inverse powers of the number of e- foldings, based on 
our experience with two-field models. 

D. Discussion 

We saw in the previous three subsections that there is little hope for a considerable amount of 
non-Gaussianity due to the evolution of modes once they cross the horizon in multi-field models 
of inflation. We restricted ourselves to separable potentials when we derive the expression for the 
nonlinear parameter in (|35l) . and looked at a few specific models with potentials of the m 2 tp 2 type: 
in the case of two fields, the contribution due to the *4-matrix is suppressed by the number of 
e-foldings, with an exponent given by twice the ratio of the heavier square mass to the lighter one. 

Extrapolating this result to the multi-field case, we concluded that a narrow spectrum would 
be the most promising candidate for a large non-Gaussian signal. This, and the fact that the most 
useful multi-field models have a narrow spectrum, eg. in the case of assisted inflation, lead us to 
evaluate the non-Gaussianity in a general multi-field scenario with a narrow mass distribution. As 
expected, we arrived at an expression that is suppressed by the number of e-foldings, in agreement 
with our experience from the two-field case. What is more, no enhancement due to the potentially 
large number of fields was found. Finally, the expression scaled like 5 2 , where 5 is a measure of the 
mass distributions width (the larger 5, the broader the spectrum). 

For the AA-field case, we used equal energy initial conditions only, for reasons of simplicity and 
since we saw in the two field model that only the numerical prefactor gets altered by a different 
choice. Consequently, we expect the general scaling behavior and the suppression by the number 
of e-foldings to be generic features, independent of the chosen initial field values. Nevertheless, 
one should perform a more careful analysis of initial conditions for any concrete model of the early 
universe, if the model is based on multiple fields with quadratic potentials (like N-flation). The 
machinery for evaluating a non-Gaussian signal developed in this paper can of course be used. 

We did not consider more intricate potentials in our case studies, like A<^ 4 potentials, or broad 
mass spectra, such as in N-flation, since our main aim here was to develop the general formalism 
and to demonstrate its applicability. We do not expect drastic new features for different potentials, 
but a thorough analysis of models in the literature should be performed, including the computation 
of the first term in (|35p 311 ] . 



It should be noted that we have restricted ourselves to slowly rolling fields throughout. However, 
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if the mass spectrum is broad, fields will leave slow roll early on, while inflation still commences. 
These fields will lead to an additional production of non-Gaussianities, due to the conversion of 
isocurvature modes to adiabatic one s 1231] ; based on the simple models studied in the literature we 

23], but more intricate models should be examined carefully 



expect these signals to be transient 
before drawing general conclusions. We also neglected any cross 



32]. 



coupling terms between fields, 



which might be yet another source of non-Gaussianity, see e.g. 

Before we conclude, we quickly derive the general expression for the scalar spectral index in case 
of a separable potential, that includes the evolution of modes once they cross the horizon. Since 
this expression incorporates the effect of isocurvature modes, it might be quite useful whenever 
these modes are present, e.g. in models incorporating compactifications of higher dimensions. 



IV. SCALAR SPECTRAL INDEX 

With our Pledge of N k a^ N u , r o m sectionHHwe caa easilv evaluate the scalar spectral 
index in the 5N formalism. In [16|] the general expression for was derived as 

" c -^- 2£+ ^ ' <98) 

n 

which can be shown to be equivalent with the result of Sasaki and Stewart [20] within the slow roll 
regime. If we insert our expressions for the derivatives of the expansion rate from (|29p and (|32p 
we arrive at 

n c - 1 = -2e* - 4 V * ± ~ > (") 

h,k=l el 

where we made use of the definition y / 2e|' = W^2>H^ip* and the Friedman equation 3H% = W*. 
It is easy to see that the last sum in the numerator vanishes, since we established already in (|4ip 
that the sum over a row or column of the ,4-matrix vanishes. Henceforth, the scalar spectral index 
simplifies to 

n ( -l = -2e* - 4 ^ j-i — '- . (100) 

EM ^ 
fe=l el 

n 

This result reduces to the one of [161 ] in the case of two fields. It is interesting to note that the 
.A-matrix, which was due to dZ^/dip^, cancels out of the expression for — 1. Obviously, this 
matrix is not constrained by observations of a nearly scale invariant spectral index. 



22 



This expression includes the effect of possible isocurvature modes and might be quite useful for 
models where such modes arise. We postpone a detailed study of (llOOp 's application to concrete 
models of the early universe to a forthcoming publication. 



In this article, we derived the general expression of the nonlinear parameter for separable po- 
tentials, within the slow roll approximation and including the evolution of perturbations once they 
cross the horizon. 

This formalism was then applied to several specific models with quadratic potentials. We find 
that non-Gaussianities are suppressed by the volume expansion rate. The power of the rate is twice 
the ratio of the heavier mass squared to the lighter one squared in two-field models. 

Extrapolating this result to multiple fields, we expect the suppression to be stronger for broad 
spectra than for narrow spectra. As a consequence, we focused on narrow spectra only for the 
general multi-field case. We recover the expected quadratic suppression with respect to the volume 
expansion rate, but we are also able to derive the explicit expression for the nonlinear parameter 
in terms of the width of the spectrum. 

Based on these case studies, we expect that multi-field models of inflation which are well de- 
scribed by the slow roll approximation cannot generate a large non-Gaussian signal. Consequently, 
we cannot use /jvz, to discriminate between multi-field models of assisted inflation and their sin- 
gle field analogs. Nevertheless, given a concrete multi-field model one should compute its non- 
Gaussianity, to verify that the general expectation is indeed true for the model at hand, since a 
dependence on the mass spectrum and, in case of quadratic potentials, the initial field values is 
present. Furthermore, we did not examine more intricate potentials, such as quartic or exponential 
ones, which clearly warrant further study [31! ]. 




As a bonus, we also computed the general expression of the scalar spectral index, again including 
the evolution of perturbations once they cross the horizon. This result should be applicable to 
computations of the effects of isocurvature modes. We postpone an application of this result to a 
forthcoming publication. 



V. CONCLUSIONS 
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